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The self-excitation of magnetic field by a spiral Couette flow between two coaxial cylinders is 
considered. We solve numerically the fully nonlinear, three-dimensional MHD equations for magnetic 
Prandtl numbers Pm (ratio of kinematic viscosity to magnetic diffusivity) between 0.14 and 10 and 
kinematic and magnetic Reynolds numbers up to about 2000. In the initial stage of exponential 
field growth (kinematic dynamo regime) , we find that the dynamo switches from one distinct regime 
to another as the radial width Srs of the magnetic field distribution becomes smaller than the 
separation of the field maximum from the flow boundary. The saturation of magnetic field growth is 
due to a reduction in the velocity shear resulting mainly from the azimuthally and axially averaged 
part of the Lorentz force, which agrees with an asymptotic result for the limit of Pm ^1. In 
the parameter regime considered, the magnetic energy decreases with kinematic Reynolds number 
as Re""'**, which is approximately as predicted by the nonlinear asymptotic theory (~ Re~^). 
However, when the velocity field is maintained by a volume force (rather than by viscous stress) the 
dependence of magnetic energy on the kinematic Reynolds number is much weaker. 
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I. INTRODUCTION 

The screw dynamo is a system where magnetic field is 
generated by the (laminar) flow of an electricaUy neutral, 
but conducting fluid with helical streamlines, i.e. 



{0,rn,u^) 



(1) 



in cylindrical polar coordinates (r,ip,z), with 17 and 
the angular and axial velocities, respectively. It is one 
of the simplest dynamo systems known and the most 
symmetric one in the sense that the flow can be steady 
and uniform in the azimuthal and axial directions. As 
first shown by Lortz |^] and Ponomarenko such 
a fiow can generate magnetic fields via dynamo action, 
i.e., without any external electromotive forces. Since the 
magnetic Reynolds number required for magnetic field 
generation by the screw dynamo is relatively low, this 
type of flow has been used in a series of laboratory dy- 
namo experiments in Riga (e.g. Refs. |^) which have 
recently achieved magnetic field growth and saturation 
1^, 0] . There are further plans to perform a dynamo ex- 
periment based on a similar (but time-dependent) flow 
j^, ^. Dynamo action of this type can also occur in 
the cooling systems of fast breeder reactors [|lO|. A re- 
lated successful dynamo experiment is the Karlsruhe liq- 
uid sodium facility which involves an ensem- 
ble of spiral flows. Since the magnetic Reynolds num- 
bers achievable in laboratory flows are never very high, 
it is important to understand quantitatively the excita- 
tion properties of the system and to predict measurable 



characteristics of the dynamo including the strength, lo- 
cation and time dependence of the magnetic field in the 
nonlinear regime. 

Other possible sites for screw dynamo action are astro- 
physical jets where a helical flow capable of dynamo 
action can arise from the axial ejection of plasma from a 
rotating accretion disc [p^ . 

A discussion of the screw dynamo in a broader context 
of slow dynamos was presented by Soward [p^ . Gilbert 
and Ponty [ p^ generalized the idea of the screw dynamo 
to certain non-axisymmetric flows and Ponty et al. 
applied this approach to hydrodynamically unstable Ek- 
man layers. 

In the present paper, we explore nonlinear states of the 
screw dynamo in the spiral Couette flow of a viscous fluid 
between two coaxial cylinders. Both the screw dynamo 
itself and the flow are simple enough to allow detailed 
analysis of the nonlinear behavior — a rare feature among 
MHD dynamo systems. In particular, this allows one 
to assess many of the empirical and heuristic arguments 
often applied to more complicated dynamo systems, such 
as the relevance of the marginally stable linear solution 
for the description of nonlinear states, and to understand 
the nonlinear states in considerable detail. 

The plan of the paper is as follows. We briefly review 
previous studies of the screw dynamo in Sect. ||, and then 
describe our model in Sect. [II. Our results are presented 



in Sect. IV A for t he (k inematic) stage of exponential 
growth and in Sect. IV B for saturated, non-linear states. 
The results are summarized in Sect. ^ 



II. THE SCREW DYNAMO 
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The kinematic behavior of the screw dynamo, includ- 
ing that in the spiral Couette flow, is well explored us- 
ing both asymptotic analysis [|8[ ^ |2l|, and 



2 



numerical modeling 
independent velocity field^ 




[2^. Consider a time- 
) where both the angular and 



axial velocity are functions of cylindrical radius alone, 
n — ^(r) and Uz = Uz{r). The evolution of the magnetic 
field B is governed by the induction equation 



— = V X (u X B - 77V X B) . 



(2) 



supplemented with appropriate boundary conditions. 
Here 77 is the magnetic diffusivity, which is related to the 
electrical conductivity a of the medium as rj = l/(/^oO')- 
At the kinematic stage, when the magnetic field is weak 
enough, u can be considered fixed and independent of B. 
The magnetic field can then grow exponentially provided 
the magnetic Reynolds number is above a certain critical 
value Rn\ \ and Eq. becomes an eigenvalue problem. 
The field is necessarily non-axisymmetric (in accordance 
with Cowling's theorem, e.g. Ref. ||2^) and, due to the 
symmetry of the flow, is a superposition of eigenmodes 
given in cylindrical polar coordinates by 

i?,(r,^,z,t) = 4-(r)e'(™^+'=^)+^*, J=r,^,z, (3) 

where m and k are the azimuthal and axial wavenumbers, 
respectively, and 

A = 7 + icj 

is the eigenvalue, with 7 the growth rate and w the os- 
cillation frequency of the magnetic field. 

While Ponomarenko Q discussed a rigid cylinder mov- 
ing in a conducting medium — thus giving rise to a dis- 
continuous velocity profile — later models |l^ |l^ 
apply to more realistic, continuous velocity fields like the 
spiral Couette-Poiseuille fiow, of which the spiral Cou- 
ette flow used in the present paper is a special case. 

The coupling of the radial and azimuthal components 
of Eq. (||), required for the magnetic field to grow (7 > 0), 
occurs via the diff usion term and is thus proportional to 
•q — see Eqs. ( |Al[) and (A2). Therefore, the growth rate 
of any given magnetic eigenmode (i.e., for fixed m and 
k) tends to zero as 77 ^ 0. The scaling of the growth 
rate with the magnetic Reynolds number cx r/^^ de- 
pends on the flow properties. In the asymptotic limit 

— 1/2 

Rm ^ 1, 7 = 0{Rm ) for a continuous velocity field 

[ p8[ , whereas 7 = 0(i?m^^^) for a discontinuous veloc- 
ity field |2^. The eigenfunction has a maximum at a 
radius ro where the advection term mfl(r) + kuz{r) [see 
Eqs. (Al) and (A2) in Appendix A] has an extremum 
in r, thus minimizing destruction of the magnetic struc- 
ture by the r-dependent advection. This implies that tq 
satisfies 



mn'iro) + ku'ziro) - 0, 



(4) 



where primes denote the derivative with respect to r. 
(In a discontinuous flow, the eigenfunction is localized 
at the discontinuity.) Modes with different ratios k/m 
are localized at different radii. An additional necessary 



condition for the existence of growing modes, which is 
due to |n9|, requires that 



d\n\n'/u'z\ 



dlni 



< 4 



(5) 



is satisfied at r = ro. This is always the case for the 
spiral Couette ffow. 

The oscillation frequency of a mode localized at r = ro 
is dominated by the advection term, 

to = -mn{ro) - fcu.(ro) + OiR;^'/^), 



for a continuous ffow with m,k ^ ^(1)- 

(crl 

The critical magnetic Reynolds number i?m , above 
which 7 > 0, depends on the radial velocity profile and is 
about 20 or larger ||, |26| . The field concentrates in 
a cylindrical shell of width Sr = 0{Rm^^^) (for R^ ^ 1) 

1 — —1/3 

for a continuous velocity field M and 6r = C'(i?m ) for 
a discontinuous velocity profile |||], provided m,k — 0(1). 
At distances from tq larger than Sr, advective distortion 
of the nonaxisymmetric magnetic field cannot be bal- 
anced by local dynamo action. Therefore the magnetic 
field must be weaker than in the resonance shell around 
ro and decays exponentially in (r— rp)^. Gilbert Q 
has obtained asymptotic solutions for the fastest mode. 



for 



m,k = 0(i?m'^) for continuous, and m,k = 0(i?m^) 
discontinuous velocity fields. 

The nonlinear behavior of the screw dynamo has been 
studied only recently in a paper by Bassom and Gilbert 
p9| , who have carried out an asymptotic analysis of the 
nonlinear case in the limit Re ^ R^n ^ 1 , where Re is the 
kinematic Reynolds number. This implies a small mag- 
netic Prandtl number, Pm = i?ni/R.e <C 1. The basic idea 
of their approach is that the overall effect of the magnetic 
field on the ffow is dominated by the azimuthally and ax- 
ially averaged Lorentz force. The exponential growth of 
the kinematic stage is saturated via a reduction in the 
velocity shear in the vicinity of r = rg where dynamo ac- 
tion is most efficient. In the asymptotic limit considered, 
the velocity shear is fully suppressed (to a given asymp- 
totic order) by magnetic forces in a shell of a radial width 

~ 0(i?n/^^"). Outside the shell, where the magnetic field 
is weaker, magnetic diffusion and stretching balance each 
other to maintain the magnetic field against Ohmic de- 
cay. At still larger distances from ro, shear dominates 
and the magnetic field is weak as at the kinematic stage. 
The steady-state field strength B^s in the spiral Couette 
ffow is estimated as IM 



0«/5Re-i), 



(6) 



where U is a characteristic value of the velocity and /xq 
denotes the magnetic vacuum permeability. As we argue 
below, the scaling with Re is sensitive to the nature of 
the driving force and arises in Eq. (^) because the ffow is 
driven by viscous forces. 
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FIG. 1: Geometrical configuration of the simulations. The 
outer cylinder rotates at angular velocity Q2, while the inner 
cylinder moves in the axial direction at speed Wi. 



Nunez et al. ijS^l have performed classical perturbation 
analysis close to the dynamo threshold for a system which 
is based on Ponomarenko's discontinuous dynamo model 
1^, but assuming the interior of the cylinder to contain 
fluid of fixed viscosity. For this sufficiently different sys- 
tem, they also find the asymptotic scaling B^^ ^ Re~^, 
while no conclusions on the asymptotic dependence on 
i?ni can be drawn from such a model. 



III. THE MODEL 

A. Spiral Couette flow 

The geometry of our model is shown in Fig. |l|. The 
conducting fluid is confined in the gap Ri < r < R2 
between two coaxial cylinders that move with axial ve- 
locities Wi and W2 and rotate with angular velocities fli 
and Q2, respectively. We choose fli = W2 = 0; the result- 
ing flow then trivially satisfles Rayleigh's stability crite- 
rion, ri2^2 > ^iRi ||3l|. The magneto-rotational insta- 
bility of the Couette flow is discussed in Refs. [^2], 
our fiow is stable with respect to this instability because 
dn/dr > 0. 

In the absence of a magnetic field, viscosity causes the 
fluid between the cylinders to adjust itself to the spiral 
Couette velocity profile 



,(C) 



where 



Ci = 



d2R'> 



-Ri 



C2 = 



C2ln:^, (7) 
r 



Wi 



The velocity profile (0), driven by the viscous stress, 
adjusts itself over the viscous relaxation time. 



{R2-R1T' 



(8) 



B. Basic equations 

The equations we solve are the MHD equations for the 
vector potential A, the velocity u and the density g: 



dA 
'di 
Du 

Dg 
Dt 



— = uxB + ?7AA + (V- A)V?7, 



--Wp+- (Au+ iVV • u) 

g g ^ ' 



jxB 



(9) 
(10) 

(11) 



complemented by an isothermal equation of state, p — 
c^g, with constant speed of sound Cg. Here, the mag- 
netic flux density B and electric current density j are 
given by B = V x A, ^oj = V x B; we denote with 
D/Dt = d/dt+ {u- V) the advective derivative; fj, is the 
dynamical viscosity (assumed constant). Below we refer 
to the Reynolds number based on the average kinematic 
viscosity, v = u/'g. 

Equation (|h implies the gauge 



77V • A + $ = 0, 



(12) 



ln(i?2/i?i 



where $ is the electrostatic potential, related to the elec- 
tric field E by E = -V<J> - dA/dt. This gauge proved 
to be most convenient for numerical purposes. Equa- 
tions (||)-(|l|) are written for compressible fiuids, but our 
choice of parameters makes compressibility insignificant 
since the speed of sound is a factor of two larger than 
the maximum fiuid velocity, which results in a density 
contrast of < 12%. 

We use an explicit finite-difference scheme of sixth or- 
der in space and third order in time described, e.g., in 
Ref. The velocity field outside the fiuid shell, i.e. 

for r < Ri and r > i?2, is prescribed and fixed, with 
u (0, 0, Wi) in r < i?i and u = (0, na?", 0) in r > R2. 
We embed the cylinders into a Cartesian box and solve 
Eqs. (|9|)-(pTl) on a Cartesian mesh in order to avoid a 
coordinate singularity on the axis and to retain the ap- 
plicability of the code to different geometries. 

The magnetic diffusivity rj is assumed to be constant 
for r < R2 — S5x (with Sx the mesh size), i.e. the inner 
cylinder has the same electric conductivity as the fluid, 
but 77 smoothly decreases to zero in R2 — 3Sx < r < R2. 
Thus, the last term in Eq. (||) is only relevant close to 
the outer boundary of the fluid. The outer cylinder is 
assumed to be magnetically impenetrable, which confines 
the magnetic field to the region r < i?2- This would 
best be achieved with a perfect conductor at r > i?2, 
but this corresponds to an infinite magnetic Reynolds 
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TABLE I: Parameters of the numerical simulations discussed in the text. For all runs, _Ri — 0.3, R2 = 1.2, = 4, Qi — 0, and 
W2 = 0. The values of Wi, 5^2 , and rj different from those of Model la are highlighted with bold face. The mesh has equal 
spacings in all three dimensions, 5x = 5y — Sz in all models except Model li where 5y = 5x = 5z/2. Lengths are measured in 
units of -Ro, angular velocities in f2o, viscosity and diffusivity in R^Q.^ and velocities in RoQ,o. 



Model 


Wi 


^2 


V 




U 




Re 






&x 


la 


1.0 


1.3 


2 X 10"^ 


7 X 10"^ 


1.85 


1.56 


111 


318 


2.9 


0.033 


lb 


1.0 


1.3 


7 X 10"^ 


7 X 10"^ 


1.85 


1.55 


31.8 


318 


10 


0.067 


Ic 


1.0 


1.3 


1 X 10"^ 


7 X 10"^ 


1.85 


1.55 


222 


318 


1.4 


0.067 


Id 


1.0 


1.3 


2.5 X 10"^ 


7 X 10"^ 


1.85 


1.55 


889 


318 


3.6 


0.033 


le 


1.0 


1.3 


1 X 10"^ 


7 X 10~^ 


1.85 


1.55 


2220 


318 


0.14 


0.017 


If 


1.0 


1.3 


2 X IQ-^ 


8 X IQ-^ 


1.85 


1.55 


111 


278 


2.5 


0.067 


Ig 


1.0 


1.3 


2 X 10^2 


4 X IQ-^ 


1.85 


1.55 


111 


555 


5.0 


0.067 


Ih 


1.0 


1.3 


2 X IQ-^ 


2 X IQ-^ 


1.85 


1.56 


111 


1110 


10 


0.033 


li 


1.0 


1.3 


2 X 10-2 


1 X IQ-^ 


1.85 


1.56 


111 


2220 


20 


0.017 


2a 


0.5 


1.3 


2 X 10^2 


4 X IQ-^ 


1.64 


1.55 


98 


491 


5.0 


0.033 


2b 


0.5 


1.3 


1 X IQ-^ 


4 X IQ-^ 


1.64 


1.55 


197 


491 


2.5 


0.067 


2c 


0.5 


1.3 


7 X IQ-^ 


4 X IQ-^ 


1.64 


1.55 


197 


491 


18 


0.067 


2d 


0.5 


1.3 


2 X IQ-^ 


2 X 10"^ 


1.64 


1.55 


98 


984 


10 


0.033 



number, the numerical implementation of which leads to 
fundamental difficulties. Therefore we use the stronger 
require ment A = for r > R2 instead. We demonstrate 
in Sect. IV A that our results are consistent with those 



obtained with a perfectly conducting outer cylinder since 
the magnetic field tends to concentrate close to the inner 
radius i?i and thus the outer boundary condition only 
weakly affects the solution. We have counter-checked our 
results with a modified magnetic condition, where the 
vector potential is 'softly' set to zero in the region r > R2 
by means of an additional relaxation term — A/ta in the 
induction equation (^), and we only report results that 
are not qualitatively affected by this change. 

In all three directions, periodic boundary conditions 
are assumed, which are imposed on the faces of the com- 
putational box. The horizontal boundary conditions are 
not actually important, since both the fluid and the mag- 
netic field are confined to r < i?2- In the axial di- 
rection (the z-direction), the assumption of periodicity 
introduces a maximum wavelength (the vertical size 
of the box) and leads to a quantization of the vertical 
wavenumber k to 



2Trn/Lz 



(13) 



with integer n, for solutions that are harmonic in z. To 
assess the role of this quantization, we have carried out 
one numerical experiment with a four times larger length 
~ 16, which introduces additional modes, which can 
have values of k four times closer to each other. The 
result was qualitatively quite similar to our other models. 
The only significant difference was that the system now 
had two modes with different fc, but very similar growth 
rates, which introduced a relatively slow drift of energy 
from one mode to the other. We will not discuss this 
simulation in the following, but rather focus on the case 



of one fixed value of to ensure the comparability of 
the different models. 

Figure ^ shows how the quantization of k due to the 
finite value of affects the dynamo system. We show 
the dependence of the kinematic growth rate 7 on the 
continuous wave number k — obtained by solving the 
eigenvalue problem ( A1)-([A^) as described in Sect. IV A 



— and indicate the quantized values kn which occur in 
our simulations. While the maximum growth rate can 
be up to 40% larger than the maximum rate measured 
at fc„ (Model Id), the optimal values of k do not differ 
much from k\. Thus, our choice of does not impose 
unrealistically small vertical scales on the magnetic field. 



C. Parameters 

The dynamo and hydrodynamical properties of the sys- 
tem are characterized by two non-dimensional quantities, 
the kinematic and magnetic Reynolds numbers, which 
are defined as 



Re = UR2IV, 



Rn 



UR2/V, 



(14) 



where U = vXR2^^2p+W^ is a characteristic velocity. 
The ratio of the two Reynolds numbers is the magnetic 
Prandtl number 



P„, = i?n./Re. 



(15) 



We introduce a reference radius Rq and a reference angu- 
lar velocity ^lo such that R2 = 1.2i?o and il.2 = f -Silo, i.e. 
the radial extent of the system and the revolution time 
are of order unity when measured in units of Rq and f2o j 
respectively. In the following, we will measure length in 
units of Rq, time in units of l/flo and velocity in units 
of RqQq, without explicitly indicating this in the text. 
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FIG. 2: Kinematic growth rate 7 versus vertical wave number 

k for Models If (• ■ • ), la ( ), Ig (- -), Ih ( ). and 

li (—■.■ — ■.■). The first two wave numbers fci, k2 allowed 
hy — 4 are indicated by vertical arrows. 



The parameter range investigated here is indicated in 
Table Parameters that remain constant for all the 
simulations are the cylinder radii = 0.3 and R2 — 
1.2), and the size of the computational box: Lz = 4 
for the vertical size and = Ly = 2(i?2 + ?>5x) for the 
horizontal dimensions. The numerical resolution is either 
42x42x60 (corresponding to 5x « 0.067) or 78x78x120 
grid points {5x w 0.033) except for Model li where the 
higher resolution of 144 x 144 x 120 {5x w 0.017) was 
used. 

For the parameters used, the viscous relaxation time 
is in the range 1-80 for the models discussed below. 
For comparison, the growth time of the magnetic field is 
greater than about 43 for the models presented, and the 
turnover time is 27r/|il2— ~ 4.8. 

Models 1 and 2 differ in the value of Wi and have 
therefore different velocity profiles. This results in dif- 
ferent dynamo efficiencies; for example, Rnl^ — 218 in 
Model 1 and i?^^''^ = 384 in Model 2. Another important 
distinction is the position of the magnetic field maximum, 
tq: in Model 1, the magnetic field is localized closer to 
the inner cylinder than in Model 2 (see e.g. the insets 
in Figs. 1^ and ^ or Fig. |l^ in Appendix B). The models 
are further subdivided (la-f, 2a-d) according to the val- 
ues of V and 77 in order to explore the effects of varying 
magnetic and kinematic Reynolds numbers. Most of our 
models have > 1, except Model Id (P,„ w 0.36) and 
Model le (P^ = 0.14). 



IV. RESULTS 

A. Kinematic regime 

The velocity profile (^ can be considered as fixed, and 
Eq. (13) as linear in B, as long as the magnetic stress is 



0.010 



0.001 







1 1 1 ' - 


- 








1 ' 






\ 1 






r 






1 


A fi„,=2-10' 














; 






7, 




; , , - 



10^ 



10-^ 10^ 



10" 




FIG. 3: Kinematic growth rate 7 and frequency u in Model 1 
for the mode m = 1, A; = fci (the dominant mode) as a func- 
tion of the magnetic Reynolds number, obtained from the 



one- dimensional eigenvalue problem described in Sect. [VA 



Dotted lines show the dependencies 7 oc Rm^^^ and 7 oc Rm'' 
in the upper panel and (u) — const) oc R^^^^ in the lower 
panel. The asymptotic solution given in Appendix^ is shown 
dashed. Results from the three-dimensional simulations are 
labeled according to the grid spacing: 5x — 0.13 (+), 0.067 
(o), 0.033 (□), and 0.017 (*). The insets show radial magnetic 
energy profiles in the domain < r < 7?2 for Jim = 2 x 10^ 
(top) and Rm = 2 x 10^ (bottom); the region occupied by the 
inner cylinder is shaded. The generation threshold for this 



configuration is R 



(cr) 



218. 



weak compared to the viscous stress, 

U 

— « ^"IT-R-- 



(16) 



Equations for the resulting kinematic dynamo problem 
are given in Appendix ^ As discussed in Sect. |l[ they 
represent a one-dimensional eigenvalue problem which is 
relatively straightforward to solve numerically. For Mod- 
els 1 and 2 we show in Figs, ^and ^the dependence of the 
growth rate 7 and frequency co on the magnetic Reynolds 
number (solid lines) and compare them to the asymptotic 
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FIG. 4: As in Fig. bl but for Model 2. The critical magnetic 



Reynolds number is Rm ~ 384 



the radial width of the field distribution decreases with 
Rm and eventually the field at r = i?i becomes negligible 
(see inset in the bottom frame of Fig. and the scaling 
7 ~ i?m is recovered. On the other hand, the field is 
always small near the boundaries in Model 2 (see insets 

in Fig. and so the scahng 7 ^ Rm^^^ does not occur. 

We have also explored the linear stage of magnetic 
field evolution using the three-dimensional code in order 
to assess its performance. A discussion of the numer- 
ical aspects is given in Appendix ^ We initialize the 
simulations with a weak random magnetic field. After 
the initial transients have died away, exponential growth 
of magnetic energy is established, corresponding to the 
fastest growing mode. 

Figure |^ shows the magnetic field structure for 
Model Ih. The level surfaces of |B| have the form of two 
helical flux tubes of opposite field orientation, which cor- 
responds to an azimuthal wavenumber m = 1. The ver- 
tical wavenumber of the solution shown is k = 2'it/Lz — 
TT /2. This mode is excited in all the models of Table |. 
However, additional higher modes are excited in Mod- 
els Ih, li, and 2d, where i?ni is larger (see Table 

Both the flow and the magnetic field are strongly heli- 
cal, the two helicities being of opposite sign (the stream- 
lines are right-handed spirals, while the magnetic field 
lines form left-handed helices as can be seen in Fig. ||). 
Since the screw dynamo mechanism relies on diffusion, 
the approximate conservation of magnetic helicity in 
highly conducting media, which leads to serious difficul- 
ties in mean-field dynamo theory ||3^, does not lead 
to any problem here. 



B. Nonlinear Simulations 



formulae, which are given in Ref. |21[| and are also repro- 
duced in Appendix ^ (dashed lines in Figs. ^ and^). The 
insets show the radial profile of magnetic energy for two 
values of 

The growth rate 7 becomes positive at R^ = i?m , 
first quickly increases with j^m, and then decreases as ex- 
pected for a slow dynamo ^ . For the continuous ve- 
locity profile (|^), the analytic theory predicts an asymp- 
totic decrease 7 c>c Rm^'^ (Appendix and this scaling 
is indeed reached for very large values of i?„i. However, 
while the growth rate for Model 2 (Fig. ^ agrees well 
with the asymptotic result (dashed) , the agreement is not 
so good for i?m < 2000 in Model 1 (Fig. |). Moreover, 
the latter model shows a transient approximate scaling 

— 1/3 

7 ~ Rn\ ■ Incidentally, this scaling is close to that for 
a flow with discontinuous radial profile |l^ . 

The difference can be explained as follows. For moder- 
ate magnetic Reynolds numbers, the field has noticeable 
strength at the boundary of the inner cylinder in the flow 
of Model 1 (see inset in the top frame of Fig. I). There- 
fore, the asymptotic theory is inapplicable as it is based 
on the assumption that the magnetic field is concentrated 
far from the boundaries r = i?i and r = i?2. However, 



After a phase of exponential growth, magnetic energy 
levels off at a certain saturation value. The corresponding 
magnetic energy density is comparable to (but smaller 
than) the kinetic energy density in the sheared flow, as 
can be seen from Table |l|, where we compare the maxi- 
mum magnetic flux density to \/^U^, and the magnetic 
energy to E^kin = ^'qU^V , where V — LzT:{RI—RI) is the 
fluid volume. The location and width of the magnetic en- 
ergy maximum are characterized by 



rs = [r) 



Stb = {{t-tbY) 



1/2 

B2 ■ 



(17) 



where {,f)j^2 = J f B'^ dV/ J B'^ dV is the magnetic- 
energy weighted volume average. These quantities are 
similar to ro and 6r of the kinematic theory. 



1. Spatial structure of the magnetic field 

Figure H shows the vertical profiles of the magnetic field 
in the saturated regime for various magnetic Reynolds 
numbers. The vertical component of the magnetic field 
is plotted as a function of z close to the radius where the 



FIG. 5: Three-dimensional representation of the magnetic field in Model Ih. (a) t — 200 (end of the kinematic stage); (b) 
t — 912 (the saturated state). Shown are surfaces where |B| is 0.65 of the maximum and magnetic field vectors. The vertical 
component of the field has opposite sign in the two flux tubes shown (mode m — 1). The azimuthal and axial variations of the 
field are harmonic at the kinematic stage (a) and fiattened in the saturated state (b) (see also Fig. 0) . 



TABLE II: Magnetic field properties in the models of Table |l[ Numbers in parentheses refer to the kinematic stage, all 
other quantities are for the saturated state. Shown are the growth rate 7 and oscillation frequency uj of the leading mode; 
the maximum saturated magnetic field strength Bmax and magnetic energy _Emag in the flow region, both normalized to the 
corresponding kinetic quantities; and the position and width of the magnetic fleld distribution in radius, rs and Stb- The 
last column lists the growing modes in the form {'m,kLz/{2n)), Lz ~ 4, ordered by decreasing growth rate. The results for 
Models Ih, li, and 2d refer to the (1, 1) mode only. Growth rates and oscillation frequencies are measured in units of Qo, 
lengths in _Ro. 



Model 


7 








-£^mag 








5rB 


{m,2k/n) 










iSkin 










la 


(0.0160) 


-2.03 


(-1.957) 


0.33 


0.015 


0.49 


(0.47) 


0.17 


(0.16) 


(1,1) 


lb 


(0.0144) 


-2.014 


(-1.957) 


0.51 


0.034 


0.48 


(0.47) 


0.17 


(0.16) 


(1,1) 


Ic 


(0.0144) 


-2.037 


(-1.957) 


0.24 


0.008 


0.48 


(0.47) 


0.17 


(0.16) 


(1,1) 


Id 


(0.0144) 


-2.045 


(-1.957) 


0.14 


0.002 


0.50 


(0.47) 


0.17 


(0.16) 


(1,1) 


le 


(0.0144) 


-2.051 


(-1.957) 


0.08 


0.0009 


0.51 


(0.47) 


0.17 


(0.17) 


(1,1) 


If 


(0.0087) 


-2.014 


(-1.976) 


0.25 


0.008 


0.48 


(0.47) 


0.17 


(0.17) 


(1,1) 


Ig 


(0.0232) 


-2.028 


(-1.935) 


0.42 


0.024 


0.47 


(0.47) 


0.15 


(0.14) 


(1,1) 


Ih 


(0.0176) 


-2.012 


(-1.91) 


0.47 


0.030 


0.47 


(0.46) 


0.13 


(0.12) 


(1,1), (1,2) 


li 


(-) 


-1.981 


(-) 


0.55 


0.035 


0.46 


(-) 


0.11 


(-) 


(1,2), (1,1), (2,3), (2,1), (2,2) 


2a 


(0.0023) 


-1.459 


(-1.471) 


0.17 


0.006 


0.63 


(0.66) 


0.18 


(0.18) 


(1,1) 


2b 


(0.0039) 


-1.459 


(-1.473) 


0.23 


0.10 


0.65 


(0.67) 


0.21 


(0.20) 


(1,1) 


2c 


(0.0039) 


-1.506 


(-1.488) 


0.48 


0.41 


0.63 


(0.66) 


0.20 


(0.20) 


(1,1) 


2d 


(0.0071) 


-1.506 


(-1.467) 


0.29 


0.01 


0.61 


(0.66) 


0.17 


(0.17) 


(1,1), (1,2) 



field concentrates. Note that the antisymmetry of the curve with respect to the middle of the box indicates that 
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only Fourier components with odd vertical wavenumbers 
are excited, which can be understood from the struc- 
ture of the nonlinear terms in Eqs. The profile 
in the saturated state is flattened in comparison to the 
(kinematic) eigenmode and this effect becomes more pro- 
nounced as i?m increases. The flattening of the maxima 
in both z and (p can also be seen in Fig. ^, where the 
level surface |B| = 0.65 |B|max is shown toward the end 
of the linear phase (Fig. ||a) and in the saturated state 
(Fig. |b). 

The strongly anhar monic profiles B{z) fo und for large 



i?„i do not occur in Bassom and Gilbert's p9| theory. 



which predicts harmonic profiles in the saturated state. 
The difference may be due to the fact that here Pm > 1, 
whereas Bassom & Gilbert assume Pm <C 1. The sensitiv- 
ity of the solution to the value of the magnetic Prandtl 
number (if this is the true reason for the difference) is 
striking. 

As can be seen in Fig. 0, the nonlinear distortion of 
the magnetic field distribution is only prominent in the 
azimuthal and axial profiles, but not much in the radial 
profile. We also note from Tabic ^ that in all the models 
the radial width brs of the magnetic energy distribution 
in the nonlinear stage is not significantly larger than that 
in the kinematic stage. In Model Ih, for example, the 
radial width of the eigenfunction is Stb ~ 0.12 during 
the linear stage and 0.13 at saturation. These numbers 
hardly differ and are very close to Sr from the asymptotic 

— 1/3 

theory, which equals Sr = R20{Rm ) — 0.11 for a dis- 



-1/4. 



0.20 



continuous velocity profile and 5r — R20{Rn 
for a smooth profile. In all our simulations, the radial 
width of the magnetic field distribution in both linear 
and nonlinear states are very close to each other and also 
close to that predicted by the asymptotic theory for the 
kinematic dynamo. This is at variance with the nonlin- 
ear asymptotics of Bassom and Gilbert |^ that predict 
the development of a 'core' region in the radial profile of 
magnetic field with a width of order 0{Rm^^^^) ~ 0.5 in 
Model li. A possible reason for this discrepancy might 
be that our models have too small values of Pm (2220 in 
Model li and 984 in Model 2d) for the asymptotic regime 
of Bassom & Gilbert to apply, even though the linear 
asymptotics are already accurate. It is more plausible, 
however, that solutions with moderate magnetic Prandtl 
number do not develop the Rm^^^^ core and their radial 
profile is quite similar to that of the kinematic eigenfunc- 
tion (see Fig. 0). 
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FIG. 6: Magnetic field structure in the saturation regime. 
Shown is, as a function of z, the vertical field Bz along a 
vertical line at r = 0.48 for Models If, la, Ig, Ih, and li; Rm 
grows monotonically along this sequence. 
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2. Nonlinear distortion of the velocity field 

Figure || shows the radial velocity profiles for the sat- 
urated phase of Model la where the scatter of the data 
points is due to their different positions in ip and z. The 
radial velocity fluctuates around zero and is dynamically 
unimportant (even more so as perturbations in the ra- 
dial velocity can be balanced by the pressure gradient). 



FIG. 7: Radial magnetic field profile in the saturated regime 
for Models li (top; Pm = 20), and le (bottom; Pm = 0.14). 
Shown is the vertically averaged magnetic energy density as 
a function of r (dots) and the corresponding profile at the 
kinematic dynamo stage (dashed). The maximum of (B^)^ 
is normalized to one; radius is measured in units of Ro- 



However, the azimuthal and axial velocities exhibit sys- 
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FIG. 8: Radial, azimuthal and vertical velocity components 
as functions of radius in the saturated stage for Model la at 
a resolution Sx = 0.033. Note the different ranges on the 
ordinate in the panels: the variation of Ur is much smaller 
than that of u^p or Uz. The shear in both uj{r) and Uz{r) 
is reduced in comparison to that in the Couette profile u^'"'' 
of Eq. (0), which is shown as a continuous gray line. The 
location of the inner cylinder is marked in gray. The scatter 
of the data points arises from their different positions in Lp and 
z. Velocities are measured in units of Ro^lo, angular velocities 
in units of Oo and radius in units of Rq. (In units of U and 
7?2, the value Ur = 0.01 corresponds to 0.0054(7, and uj — 1.2 
corresponds to O.78U/R2, for example.) 



tematic deviations from their Couette profiles so that the 
velocity shear is reduced in the region where the magnetic 
field concentrates. It is especially clear in the case of Uz 
(Fig. H^) that the spatial scatter is smaller than the mean 
variation, so the distortion of the velocity field is axially 
symmetric and independent of z to a first approximation. 
To justify this we demonstrate in Fig. |l^ that the radial 
profiles of the averaged Lorentz force are in close corre- 



FIG. 9: As Fig. H, but for Model li. 



spondence with the deviations of the respective velocity 
components of Fig. ^ from the Couette profile. This is 
compatible with the scenario of Bassom and Gilbert , 
where, in the limit of infinite kinematic and magnetic 
Reynolds numbers, the saturation of dynamo action is 
mainly due to the vertically and azimuthally averaged 
part of the Lorentz force. As the magnetic Reynolds 
number is increased, the relative distortion of the veloc- 
ity field, and in particular the reduction of velocity shear, 
are getting more pronounced; this can be seen in Fig. |^, 
where we show the same profiles as in Fig. ^ but for a 
seven times larger value of i?ni- 

The radial width of the region where the velocity field 
is distorted away from the original Couette profile is 
much larger than that of the Lorentz force (see Fig. 
below). This happens because the flow adjusts itself to 
two separate Couette profiles at both ends of the radial 
range where the Lorentz force has distorted it, and so a 
localized magnetic field does affect the flow throughout 
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FIG. 10; Vertically averaged components of the Lorentz ac- 
celeration Fl = j X B/£) as a function of radius for Model la. 
Note that the effects of magnetic torques on the inner cylin- 
der (shaded area in Fig. lid) are ignored in our model and 
its velocity remains fixed. Velocities are measured in units of 
RqQ,o, angular velocities in units of Qq and radius in units of 
Ro- Fl is measured in units of Ro^oj radius in units of Ro- 
(In units of U and R2, the value -pLor ~ 0.05 corresponds to 
0.018J7V^2, for example.) 



the domain. It can be expected that a flow profile driven 
by a volume force will be distorted less in regions where 
magnetic field is weak; our results confirm this expecta- 
tion. 

The screw dynamo can be interpreted in terms of the 
mean- field al7-dynamo jlj, |l^, the f2-term being as 
usual the shear term rBr-dfl/dr [see Eq. (A2)]. The a- 
efi'ect is identified with a part of the diffusion term, and 



the corresponding term in the induction equation has the 
form 



d(p ^ 



(18) 



[see also Eq. ( |Al[ )]. Since the differential operator act- 
ing on is obviously not infiuenced by the magnetic 
field strength, the a-effcct cannot be affected by the 
growing magnetic field and saturation is fully due to Q- 
quenching. This is opposite to one of the standard sce- 
narios for mean-field dynamos (a-quenching), where the 
magnetic field has little infiuence on the angular velocity 
and saturation is considered to be caused by the partial 
suppression of the a-effect by magnetic fields. 
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FIG. 11: Magnetic energy in the saturated state versus the 
kinematic and magnetic Reynolds numbers for Models la-i. 
Top: -Emag/Skin as a function of Re for Models lb, la, Ic, Id, 
and le (asterisks) . The c rosses refer to a model with volume 



forcing (see Sect. [VB3). Bottom: -Emag/Skin as a function 



of -Rm for Models If, la, Ig, Ih, and li. The dashed lines show 
the asymptotic power laws given in Eq. (^); the dotted line is 
a least-square power-law fit and corresponds to a dependence 



iSmag oc Re 



3. Scaling of magnetic energy with Re and Rm 

Figure |ll| shows the magnetic energy in the satu- 
rated state as a function of the kinematic and magnetic 
Reynolds numbers. We also show the asymptotic scal- 
ings of Bassom and Gilbert [^. There is clearly some 
agreement between the numerical and asymptotic results 
although the asymptotic solution has been obtained for 
Pni 1 whereas the numerical results refer mostly to the 
case Pm > 1- We believe that the reason for the rough 
agreement is related to the fact that, in both cases, the 
magnetic feedback affects mainly the axisymmetric flow 
profile, as was seen in Fig. ^. 

As expected (at least for weakly supercritical solu- 
tions) the steady-state magnetic field strength B^s in- 
creases with magnetic Reynolds number. On the other 
hand, the saturated field strength Bgs decreases with Re 
in a fiow driven by viscous stresses because it becomes 
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easier for the magnetic field to modify the velocity field 
in a given cylindrical shell as Re increases (owing to the 
weaker viscous coupling of fluid at different radii), and 
so a weaker magnetic field is needed to achieve a local 
reduction in shear sufhcient to halt the field growth. As 
we discuss below, this behavior is characteristic of a flow 
driven by viscous stresses rather than by a volume force. 

The scaling of magnetic energy with Re is slightly shal- 
lower than predicted by Bassom & Gilbert and has an 
exponent close to —0.84 rather than —1. Similarly, the 

2 /5 

growth of magnetic energy with is slower than i?m 
at the values of i?m explored here. The most plausible 
reason for these disagreements is the difference in mag- 
netic Prandtl numbers and, possibly, also the fact that 
the magnetic field distribution is still not narrow com- 
pared to the gap width R2—R1 for the magnetic Reynolds 
numbers we were able to consider. 

Taken at face value, the dependence of Bgs on Re dis- 
cussed above implies that the resulting magnetic field 
will be negligible wherever Re ^ 1. However, in most 
real systems the fiow will be driven by non- viscous forces. 
These are pressure and inertial forces in the dynamo ex- 
periments in Riga and Perm. In the case of astrophysi- 
cal jets, acceleration and collimation of the flow can be 
due to an external magnetic field (axisymmetric to the 
first approximation, then the screw dynamo can generate 
an additional nonaxisymmetric magnetic field). In these 
cases the kinematic Reynolds number will not play such 
a prominent role in the system and Bgs is expected to be 
independent of Re for Re ^ 1 . 

To verify this idea, we have carried out further numer- 
ical simulations for the parameters of Models lb, la, Ic, 
and Id, but with an additional volume force 



u— u 



(C) 



(19) 



tually achieves energy equipartition with the turbulence. 
For small magnetic Prandtl number, however, the dy- 
namo is only weakly supercritical and only large scales 
can be excited. Thus, the screw dynamo effect is ex- 
pected to work, but it is now controlled by the effective 
values of Re and i?ni based on turbulent diffusivities. We 
conjecture that this could still be qualitatively true in 
the nonlinear phase of field evolution. Since the turbu- 
lent values of Re and R^ can be quite moderate, even 
the scaling (||) would not lead to physically uninterest- 
ing magnetic fields. It can be expected, however, that the 
turbulence in the Couette flow — being generated mainly 
in turbulent boundary layers near the flow boundaries — 
will be very inhomogeneous and this can affect the the- 
ory discussed here. On the other hand, the reduction 
in the effective magnetic Reynolds number due to turbu- 

(cr) 

lence has to be only moderate since one needs R^ > Rin 
for any sort of dynamo action. Assuming that turbulent 
magnetic and kinetic Reynolds numbers have similar or- 
ders of magnitude. Re ~ Rm, and that R^i^ ~ lO'^, the 
resulting magnetic energy density will be about 10% of 
the kinetic energy density even with the scaling (^. 

We note in this connection that the maximum ra- 
dial magnetic field component measured in the Riga dy- 
namo experiment ^ is about 6 mT, while the equiparti- 
tion field strength would be w 0.5 T. B,^ and B^ be- 
ing about five times larger than _B,. in the kinematic 
case [F. Stefani, private communication], one gets a ratio 
£^mag/£^kin ~ 5 X 10"'^. This is two ordcrs of magnitude 
larger than i?m^Re"^ ^ 10~^ and clearly implies that 
the scaling (^, which is based on the assumption of a 
viscously driven, laminar flow, cannot be directly applied 
to flows in laboratory experiments. 



on the right-hand side of the equation of motion (10). 
Here u^*-^^ denotes the spiral Couette profile (Q) and r 
is the time scale over which the flow adjusts itself to 
the Couette proflle. With our choice t = 1, we have 
T ^ Tvisc , where the viscous time Tvisc varies in the range 
1-80, cf. Eq. (^). Thus, the Couette flow proflle is now 
maintained on a dynamical time scale rather than by 
viscosity if Re 30. 

We show in Fig. ^ by crosses the resulting dependence 
of the steady-state magnetic energy on Re; the depen- 
dence is clearly much weaker than in the case of viscous 
driving (a power-law fit to the points shown has an ex- 
ponent of about —0.41). In the inviscid limit Re — > oo 
the magnetic energy must eventually become indepen- 
dent of viscosity, and so we anticipate that the depen- 
dence marked with crosses has a horizontal asymptote. 

Another reason why the dependence (|^) is not directly 
applicable to laboratory and astrophysical dynamos is 
that the corresponding flows are turbulent. In that case 
the generation of magnetic field on the scale of the mean 
flow can be accompanied by the so-called 'fluctuation dy- 
namo' producing small-scale magnetic flelds, which even- 



All our experiments were carried out for angular veloc- 
ity Q increasing outwards and, for the boundary condi- 
tions and parameters considered, the flow was found to 
be hydrodynamically stable in the sense that the Lorentz 
force led to relatively small deviations from the unper- 
turbed Couette profile (0). An interesting question is 
that of dynamo action in a spiral Couette flow unstable 
with respect to the magneto-rotational instability, i.e. for 
angular velocity decreasing outwards. Flows of this type 
— although for vanishing vertical velocities of the two 
cylinders — are currently discussed in connection with 
experimental investigations of the magneto-rotational in- 
stability If the fiow is hydrodynamically 
stable (angular momentum increasing outwards), a cer- 
tain magnetic field is necessary to obtain regular Taylor 
columns; however, these columns involve a helical veloc- 
ity field and are themselves capable of dynamo action 
Q. Thus, the magnetic field can play the role of a 
'catalyst', destabilizing the system and at the same time 
regenerating it by the resulting instability; this mecha- 
nism has been demonstrated to occur in accretion discs 
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V. CONCLUSIONS 

The numerical calculations performed here have re- 
vealed a new intermediate asymptotic regime in the kine- 
matic screw dynamo where, in a certain range of i?rn, the 
magnetic field eigenfunction is large enough at the flow 
boundary so that the standard asymptotic solutions are 
inapplicable and the growth rate of the magnetic field 

— 1/3 

scales as 7 = 0{Rm )■ This dependence is typical of a 
discontinuous velocity profile. The standard asymptotic 
scaling, 7 = 0{Rm ), is only recovered at significantly 
larger values of i^m where the eigenfunction becomes nar- 
row enough in radius as to satisfy the assumptions of the 
asymptotic theory. 

We have confirmed the result of Bassom and Gilbert 
p9| that saturation of screw dynamo action occurs via a 
reduction in the velocity shear produced by the axisym- 
metric part of the Lorentz force. We have shown that 
this also applies to the case of large magnetic Prandtl 
numbers Pm. However, the radial profile of the nonlinear 
solution is very similar to the marginally stable eigen- 
mode for the nonlinearly modified velocity field. This is 
different from the nonlinear asymptotics of Bassom and 
Gilbert, which predicts a plateau in the radial depen- 
dence. It is possible that such a plateau can only occur 
at values of the magnetic Prandtl number much smaller 
than what we have been able to simulate in the present 
work, but we have not detected any tendency towards its 
development at = 0.14, the smallest value explored 
here. 

We have demonstrated that the scaling of the steady- 
state magnetic field with kinematic Reynolds number is 
sensitive to the nature of the forces driving the fiow. 



APPENDIX A: ONE-DIMENSIONAL 
KINEMATIC DYNAMO PROBLEM 

In cylindrical coordinates {r,ip,z), the kinematic 
dynamo problem for given steady velocity field 
(0,ril{r),Uz{r)) , is given by the following non- 
dimensionalized set of equations for the two amplitudes 
Br, % [cf. Eq. (I)]: 

^ ^ 1 f . ^ 2ir7? 1 
XBr + i{mn+kuz)Br = <^ DB^ —B^ ^ , (Al) 



dn 

dr 
1 

R„ 



2im : 



DB^ + ^Br y{A2) 



obtained from the solenoidality condition, 
1 d 



dr 



{rBr 



-B^ 



ikBz = 0. 



(A4) 



Second-order asymptotic analysis (for R^ 00) for 
the spiral Couette flow has been presented in Refs. 
|]l8| , [21II . It is convenient to deflne the magnetic Reynolds 
number as 



Rm — 



RiRi \m ^ ^EiriP 



V 



(A5) 



here, which is different from the definition ( |l4[) used in 
the body of this paper. The radius ro, where the mag- 
netic field concentrates, is given by 



'lni?2- 


lni?i 


i?9 - 





ro — RiRi^ 

and the field can only be growing if 
mAr2 

<0, 



— A6) 



(A7) 



/cAiy 

where Ar2 = Oj-f^i and AI^ = 14^2-14^1- 

— —1/2 

To second order in the small parameter R^a , the 
eigenvalue A of the fastest growing mode is given (in di- 
mensional units) by 

A + i[TOil(ro) + kuz{ra)\ 
V 



= ^^(V2-l)v/|m|(l+aiiX/2 



iI-V2- 
36 



; 2 2 

k rn 



o 



(A8) 



where cti — sgn(fcAII^) denotes the sign of k/SW . As 
T] oc R'^ , the leading-order term on the right-hand side 
is ©(-Rm^^^) when measured in the units \J jRi. 

To the leading order in the magnetic field ampli- 
tudes B are given by 

"(r-ro)2" 



Bip 

Br 

B, 



exp 
+ 



-(l+crii)y |m|i?„ 



2rl 



0-2(l+crii)\/ ^B^Rr^ 



1/2 



k 



B^ 



0(R7y' 



(A9) 
(AlO) 
(All) 



where (T2 = sgnAi7. 



where 

t) - ^ ( '^^ ^^^'^ k^ (A3) 
r dr \ dr J 

is a self-adjoint, Laplace-type differential operator, and 
A is the eigenvalue. The vertical component B^ can be 



APPENDIX B: ACCURACY OF THE 
NUMERICAL SCHEME 

We assess the accuracy of our three-dimensional sim- 
ulations and the implementation of the boundary condi- 
tions by comparing solutions obtained with the 3D code 
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FIG. 12: Radial field structure in the kinematic regime as 
obtained from a three-dimensional simulation at a resolu- 
tion 5x = 0.033 (dots) and from solving the one-dimensional 
Eq. with _Ri = 0.28 (solid), (a) Model la, (b) Model 2a. 
The maximum of (B^)^ is normalized to one; radius is mea- 
sured in units of Rq. 



in the kinematic regime with those from the correspond- 
ing ID eigenvalue problem. 

The radial dependence of the (vertically averaged) 
magnetic energy density is shown in Fig. |l2| for the high- 
resolution runs of Models la and 2a at a time when the 
exponential growth has well established itself. Magnetic 
energy concentrates in a cylindrical shell of radius 
and of half- width Srs, defined in Eq. (p7[). For Model la 
(Fig. |l2|a), we find vb = 0.48 and 6rB = 0.16 (see also 
Tableau) . For comparison, we have overplotted the pro- 
file obtained from the one-dimensional model. 

Figure ^ shows good agreement between the three- 
dimensional and one-dimensional simulations everywhere 
except close to the outer boundary. For r w i?2, 
the magnetic energy in the three-dimensional simulation 
smoothly turns to zero, while in the one-dimensional cal- 
culation the tangential components of B remain signifi- 
cant for Model 2a. This is due to the different boundary 
conditions used (A = vs. perfectly conducting) and the 
close agreement of growth rates and eigenfunctions gives 
us reason to believe that these are not influenced by this 
localized deviation. 



TABLE III: Accuracy of the three-dimensional code, illus- 
trated for Model la using various resolutions Sx. The effec- 
tive inner radius R^^^ is defined in the text. Growth rate, 
frequency, localization radius and width of the magnetic field 
distribution are given for Model la at different resolutions 
(top half). The lower half of the table shows the acc urat e 
values obtained by solving the one-dimensional problem ( Al ) , 
( A2) for three different values of Ri. Lengths are measured in 
units of _Ro, growth rates and oscillation frequencies in units 
of fio- 



Model 




7 


UJ 


tb 


SrB 


la ((5x=0.067) 


0.291 


0.0144 


-1.969 


0.467 


0.162 


la (5x^0.033) 


0.301 


0.0160 


-1.957 


0.478 


0.160 


la (52;=0.017) 


0.300 


0.0162 


-1.952 


0.483 


0.161 


theory (_Ri=0.28) 




0.0140 


-1.971 


0.468 


0.158 


theory (i?i=0.29) 




0.0154 


-1.961 


0.478 


0.159 


theory (i?i=0.30) 




0.0167 


-1.951 


0.488 


0.160 



0.40 



0.30 



^ 0.20 



0.10 



0.00 
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FIG. 13: Representation of the inner cylinder (shaded) on a 
Cartesian mesh for &x = 0.067 (left) and &x = 0.033 (right). 
Overlayed are the circles of effective radius ¥y\^^ , and arrows 
representing the velocity field. All lengths are measured in 
units of i?o. 



In Table III , we compare the eigenvalues A and the spa- 
tial parameters tb and bvB of the magnetic field for the 
three-dimensional simulations with those from the one- 
dimensional model for different numerical resolutions. 
While the high-resolution simulation has only an error 
of about 5% in the growth rate 7, and 0.3% in the fre- 
quency w, the lower-resolution run yields errors of 18% 
and 0.9 %, respectively. The main source of inaccuracy is 
the angular discretization of the inner cylinder boundary 
as illustrated in Fig. If we define an effective inner 
radius as the radius of a circle enclosing the same 

area as the shaded cells in Fig. |l|, then « 0.29 

for the lower-resolution run. For an inner radius of 0.29, 
the one-dimensional model yields a growth rate that is 
considerably lower and closer to what is observed in the 
lower-resolution run (see Table |l|). This adjustment of 

(cfl) 



well 



results in better agreement of a;, tb and bvB as 
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